Excitonic pairing between nodal fermions 
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We study excitonic pairing in nodal fermion systems characterized by a vanishing quasiparticle 
density of states at the pointlike Fermi surface and a concomitant lack of screening for long-range 
interactions. By solving the gap equation for the excitonic order parameter, we obtain a critical 
value of the interaction strength for a variety of power-law interactions and densities of states. We 
compute the free energy and analyze possible phase transitions, thus shedding further light on the 
uimsual pairing properties of this peculiar class of strongly correlated systems. 



Since its foundation some fifty years ago, Fermi-liquid 
theory has been providing a convenient framework for the 
analysis of conventional metals. Nowadays, however, the 
quest for departures from classic Fermi-liquid behavior 
has become one of the central themes in modern con- 
densed matter theory. 

In a generic one-dimensional (ID) many-fcrmion sys- 
tem, an arbitrarily weak short-range repulsive interaction 
is known to completely destroy the Fermi liquid, thereby 
giving rise to so-called Luttinger behavior. In higher di- 
mensions (2D and 3D), the Fermi liquid is believed to 
be more robust, although it is not expected to remain 
absolutely stable. While in the case of short-range re- 
pulsive interactions any deviations from the Fermi liquid 
are likely to be limited to the regime of infinitely strong 
couplings, the long-ranged interactions are considered to 
be generally capable of destroying the Fermi liquid even 
at finite couplings. 

In that regard, a particularly promising case appears to 
be the class of systems which possess isolated ( "nodal" ) 
Fermi points where the quasiparticle density of states 
(DOS) vanishes. As a result, any long-ranged interac- 
tion (whether the bare Coulomb one or that mediated 
by soft collective modes proliferating near an emergent 
instability) remains largely unscreened. 

Examples of the nodal fermion systems are provided 
by degenerate semimetals, such as graphite or bismuth, 
where the conduction and valence bands touch each other 
at isolated points in the Brillouin zone. In a monolayer 
of graphite, for example, the two-dimensional fermion 
excitations with the momenta close to the nodal points 
have a linear dispersion, and, therefore, can be described 
in terms of a Dirac-like effective theoryii 

Amongst other potential candidates to the nodal liquid 
scenario are p (such as Ilea) and d (such as cuprates) su- 
perfluids where the quasiparticle gap has isolated nodes. 
Yet another example is provided by charge density wave 
(CDW) insulators, such as dichalcogenides, which are 
thought to have /-wave CDW ordern^ 

On general grounds, one might expect behavior man- 
ifested by systems supporting nodal quasiparticles to be 
markedly different from that of fermions with an ex- 
tended Fermi surface from which a bare long-range in- 
teraction would undergo strong screening. In particular, 
it is conceivable that pairing instabilities which, depend- 
ing on the sign of the interaction, occur either in the 



particle-particle (Cooper) or particle-hole (Peierls) chan- 
nels, develop differently from the conventional BCS sce- 
nario. 

Barring a handful of counter-examples, most of the 
previous studies of fermion pairing were limited to the 
simplest BCS case where the momentum dependence of 
the corresponding order parameter (except, possibly, for 
its angular dependence) would be routinely neglected. 
However, while being generally justified in the case of 
normal metals with their intrinsic screening, the BCS 
solution may no longer be applicable in the nodal liquids. 

In the present paper, we address the general problem 
of pairing between the nodal fermions with a generic dis- 
persion (X . The latter corresponds to the power- 
law DOS i'{e) cx which vanishes at the nodal Fermi 
point located at p = and is governed by the exponent 

Although below we focus on excitonic pairing pertinent 
to the case of repulsive interactions, much of the follow- 
ing discussion can be readily adapted to the case of the 
Cooper instability caused by attraction. 

In the particle-hole representation, the inverse fermion 
Green function reads 



G ^{uj, p) = cj - ^ - fs^p + TiAj, 



(1) 



where is the triplet of the Pauli matrices acting in the 
particle-hole space, n is the chemical potential, and Ap 
is the gap function which describes the anticipated spa- 
tially uniform (see the conclusions for a discussion of this 
assumption) pairing in the Peierls channel caused by a 
long-ranged repulsive potential with the Fourier trans- 
form t/(°)(q) = 50/9''" ■ 

When the nodal fermion polarization x(cj, q) is taken 
into account, the renormalized effective interaction as- 
sumes the form 



C/(^,q) 



.90 



9''" +gx('^,q) 



(2) 



At low (compared to a characteristic scale A which is of 
order the maximum span of the Brillouin zone) momenta 
and still lower frequencies {uj/q^ and q ^ 0) the one- 
loop polarization behaves as x(0, q) oc q^^^. Thus, un- 
der the condition a < P, the effective interaction (2) re- 
tains its bare form (except for a finite renormalization of 

the couphng constant go go [1 + .9ox(0, q)/9''"]^^ Ig-^o 
and an additional frequency dependence for a = (3). By 
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contrast, ior a > (3 the overscreened effective interaction 
(2) acquires a universal form governed by the fermion 
polarization, thereby giving rise to the behavior which is 
qualitatively similar to that occurring at a = 

The gap equation is given by the fi -component of the 
Dyson-Schwinger equation for the fermion propagator (1) 



(3) 



where Go(w, p) is the bare propagator with Ap = 0. 

In what follows, we focus on the regime a < (3 where, 
conceivably, the polarization operator may still become 
important at the momenta of order the upper cutoff set 
by A. We find, however, that a rapid decrease of Ap 
at large momenta (see below) makes the gap function 
robust against the inclusion of such high-energy polar- 
ization effects. By the same token, a polarization-related 
frequency dependence (retardation) of the effective inter- 
action does not have a significant impact on the solutions 
of Eq.(l) and, therefore, will be neglected from now on. 

Taking into account the above arguments, we perform 
the frequency integration in Eq. (3) and take its trace with 
fi, thereby arriving at the integral equation for the mo- 
mentum (albeit not the frequency) dependent gap func- 
tion 



Ap = EC/(0,p-q) 



^^^^(a*tK), 

2E„ ^' 



(4) 



where = y -f A^ is the quasiparticle energy, and 
the sum is taken over all the occupied states in both the 
valence [0, = —Eq) and conduction (S7 = Eq) bands. 

Notably, for any a > 0, a nontrivial momentum depen- 
dence of the integral kernel in Eq.(4) rules out the pos- 
sibility of the standard BCS solution Ap — const. While 
resisting an exact analytical solution for general values of 
a and /3, Eq.(4) can still be amenable to numerical anal- 
ysis. However, we find it more physically transparent 
to resort to an approximate analytical approach which 
generalizes that of Refs. 3 and 4 and allows one to con- 
vert an integral equation into a simpler linear differential 
one. This method elaborates on the earlier analyses^ 
of the phenomenon of chiral symmetry breaking in 3D 
relativistic fermion systems that can be thought of as a 
Lorentz-invariant analog of excitonic pairing. 

We start off by observing that, at large momenta, 
Eq.(4) can be simplified by replacing the kernel U{p — q) 
with either U (p) or U{q), depending on the relative mag- 
nitude of the momenta q and p. Next, by changing 
the integration variable from the momentum p to the 
quasiparticle energy e = Ep and repeatedly differentiat- 
ing both sides of Eq.(4), we obtain the linear differential 
equation 



d^A{e) a + 1 dA{e) aA{e) 



de2 



de 
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0, 



(5) 



where g = go/[2°-^T:°/^V{D /2)]. 

In order for Eq. (5) to be consistent with the original in- 
tegral Eq.(4), it has be complemented with the boundary 
conditions 



dA(e) 



and 



de 



A(e) + 



= 



e— inax[A,/i] 



e dA(e) 



a dE 



0. 



(6) 



(7) 



In what follows, we refer to Eqs.(6) and (7), respectively, 
as the infrared (IR) and ultraviolet (UV) boundary con- 
ditions. 

Although the above approximations can be fully justi- 
fied only at large momenta, a comparison between the nu- 
merical analyses of the original gap equation (4) and the 
relatively simple approximate solutions of the linearized 
equation (5) (see below) shows that the latter provide a 
faithful representation of the solutions to Eq. (4) down to 
the lowest energies e ~ A. 

We first consider the case of half filling (/x = 0). 
The linearized Eq.(5) can be interpreted as the radial 
Schrodinger equation of the zero energy state in the fic- 
titious (a -|- 2)-dimensional spherically symmetrical po- 
tential V{R) = gal R^^"^^ . This analogy suggests that 
for < /? — a < 1 the zero energy solution occurs for 
couplings g in excess of a finite threshold value gc- 

This prediction is supported by a direct analysis of 
Eq.(5) which, for a > /3/3, appears to have an approxi- 
mate solution of the form 



A(e) 



^^-Ay^-) sin[a>(6) + ^] 
e J sin (5 



(8) 



where the prefactor is chosen to satisfy the natural nor- 
malization condition A(e)|e=A — A, and the phase of the 
trigonometric function is given by the expression 



^ iv«5e^'--«W/4[l - (A/e)n 



(9) 



Here r = (/3 — a)/2, the phase shift S = 
tan~i[^4agA2'- - K(A)/a] for A > [K(A)/4ag]^/^'', 
otherwise 5 = 0, is determined by the IR boundary con- 
dition, and the slowly varying functions 7(e) and K{e) 
attain the following asymptotic values: 



7(A) = a/2, 7(A) = (a-f/3)/4 (10) 



and 



K{A)^a^-{a-p)'^, k(A) 



-ia-P)\ (11) 



The maximum (zero-momentum) value of the gap func- 
tion A = A(0) can be determined from the UV boundary 
condition (7), which now assumes the form 



tan[$(A)] 



4 ag-K(A)(A/A)74 
3a- P ^ag - k(A)/4 



= 0, (12) 
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where g = gA'^^ . 

For < a < /3, Eq.(12) yields 

A ^ A{~g - ~g,y/^ (13) 

for couplings greater than the critical value 

gc^-^iia + pf+TT^iP-ar] (14) 
loa 

In the marginal case a = /3 > and for g > gc — a/A the 
solution (8) takes the form 



A(e) 



a/2 sin y/ag ~ Q!2/41n(e/A) + S 
e I sin J 



(15) 

where S = tan"-'^ y^4:{g/a) — 1. 

In the case a — /3 = 1, Eq.(15) reproduces the solution 
that was obtained in Rcfs. 3 and 4. In this case, the 
gap given by Eq.(12) shows a behavior reminiscent of 
that at the Kosterlitz-Thouless (KT) transition in XY- 
symmetrical 2D systems 



Aexp 



27r-4(5 
^4a(g - gc) 



(16) 



A similar behavior has been previously found in the con- 
text of 3D relativistic chiral symmetry breaking, in which 
it was conjectured to signal the onset of a conformally 
invariant critical regime »^ The present discussion shows 
that the KT-like behavior (16) is not specific to Lorentz- 
invariant systems. Rather, it appears to be characteris- 
tic of the scale invariance of the underlying "radial" gap 
equation (5). 

Next, we introduce a finite chemical potential /i, which, 
depending on the sign, represents a finite density of either 
particles or holes. The occupation factor 0{fi^Eq) in the 
right-hand side of Eq. (4) sets a lower limit of momentum 
integration. 

A straightforward analysis shows that at e the 
solution A(e) of the gap equation levels off and becomes 
independent of energy. For A < /i we then impose the 
normalization condition A(e)|(;=^ — A, thereby arriving 
at the approximate formula for the counterpart of Eq.(8) 



A(e) 



7(e) sin [$(e) + S] 
sin [$(^) + S] ' 



(17) 



whereas for A > ^ the solution (8) remains essentially 
intact. In particular, the marginal (scale-invariant) case 
a = /3 is now described by the gap function 



A(.)^A(^) 



^N^"/2 sm{^y ag ~ a^/Alne/A + 5) 
sin{yjag - Q;2/41n^/A + 6) 



(18) 



After having found the approximate solutions to the gap 
equation, we now proceed to elucidate the nature of the 
quantum phase transition which results in the opening of 



the gap. To that end, we use Eqs.(8) and (18) to compute 
the free energy functional i^(A,/i). 

The presence of the long-ranged forces makes it rather 
challenging to utilize the conventional Luttinger-Ward 
approach (for its recent application to a very detailed 
analysis of the non-BCS energy-dependent superconduct- 
ing pairing, see Ref. 6). Specifically, a straightforward 
implementation of this method gives rise to a strongly di- 
vergent double momentum integral X)p q ApC^~^(0, p — 
q)Aq whose kernel is given by the (generally, highly non- 
local) inverse operator f7~^(0,p). 

Moreover, an unequivocal identification of the solu- 
tions to the gap equation as minima (as opposed to 
maxima) of i^(A, jS) requires one to compute this func- 
tional for all A, whereas the previous applications of the 
Luttinger-Ward and similar techniques have been limited 
to the evaluation of the free energy f (Aq, /i) at its abso- 
lute minimum Aq ("condensation energy").® 

As an alternate approach, we employ a field-theoretical 
technique which has been previously used to investi- 
gate the conjectured excitonic instability in pyrolytic 
graphite.'* This method makes use of the source function 
(see Ref. 4 and references therein) 



J(A) 

and the order parameter 
ct(A) = Tr[fi 



a ae 



q 



e 

-A{e) 
90 



(19) 



(20) 



which both determine the sought-after (regularized) ef- 
fective potential given by the formula 

r-A 



F(A,/.)= / dA 
Jo 



, da{A'] 
dA' 



J(A')+F(0,/i). (21) 



The integration constant -^(0, ^) is to be found from the 
equation for the density of excess particles 

dF{A,fi) 



n{ti) 



0{pTEp). 



(22) 



In order to illustrate this approach on a well-known ex- 
ample, we first consider the standard BCS pairing be- 
tween fermions with a short-range interaction potential 
and a finite density of states {a = (3 = 0). 

For fi — and an arbitrarily weak repulsion, the gap 
equation (4) yields the conventional BCS-type solution 
A(e) = Ao = Ae-i/9. 

For /i < A the free energy reads 

whereas for A < /i one obtains 



, . , oA^ I 1 

igo \g 2 



A 



250 



(24) 
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the last term being proportional to the particle density 
n = iv'M'-A2. 

The free energy given by Eqs. (23,24) and plotted in 
Fig.l exhibits a minimum at A = Aq for all ^ < fXc ^ 
and g < 2. Upon increasing /i past the critical 
value ^c, the minimum jumps discontinuously to the triv- 
ial solution A = 0, consistent with the anticipated first 
order nature of the transition in question, which is known 
to be equivalent to the superconducting BCS transition 
in an external magnetic field. 

In the case of Fermi liquids with screened repulsive 
interactions, the same scenario describes the excitonic 
instability. It should be noted, however, that the transi- 
tion driven by a variable density (which provides a more 
appropriate picture of the system, in which the number 
of carriers is not fixed) can still be of second order. ^ 

This given example demonstrates that by merely eval- 
uating the free energy at the extremal points of F{A, /i), 
one may not always be able to distinguish the true abso- 
lute minimum from a local one or, worse yet, from a max- 
imum. In that regard, it is worth mentioning the appar- 
ent confusion over the nature of the so-called "nontrivial" 
solution A = ■\/2Ao/i — Aq that was claimed in Ref. 7 
to exist for /i > Ao/2. This solution was also invoked 
in some of the recent studies of the conjectured (weakly 
ferromagnetic) excitonic insulating state in hexaborides.* 

However, it can be immediately seen from Fig.l that 
the extremal point in question is, in fact, a local max- 
imum of the BCS free energy functional. Therefore, it 
could not have provided a viable solution to the problem 
of excitonic pairing in a doped semimetal even if the hex- 
aborides belonged to this class of materials (which is now 
considered rather unlikely from the experimental point of 
view^) . 

As the next example, we discuss the case of a short- 
range potential and a linear density of states (a = 0, 
/? = 1). Similar to the BCS case, the constancy of the gap 
function makes the solution to the gap equation [A(e) = 
Aq = g^^ — g^^ for /i < A] trivial. 

For /i < A we obtain 



9 

250' 



(25) 



while for A < /i the free energy incorporating an excess 
electron density n(/i) — ^(/i^ — A^) takes the form 



250 



1 



2^ 2 " 



;aoA2 



(26) 



Notably, the terms of order A^A which are present in 
both the kinetic and potential contributions cancel each 
other, and the condensation energy is governed by the 
fermion states with energies e ~ A. 

As follows from Eqs. (25,26), the minimum of the free 
energy plotted in Fig. 2 exists for all g > gc = 1/A. Inter- 
estingly, at fif, = Aq the mean-field free energy becomes 
independent of A in the entire range < A < Aq, thus 
signaling the possibility of a continuous transition. It is. 



however, conceivable that upon accounting for the fluc- 
tuations about the mean-field solution one would find 
the transition in question to be of (possibly, weakly) first 
order. 

Being primarily concerned with the properties of the 
systems governed by long-ranged interactions, we now 
move to the case of a linear density of states and the un- 
screened Coulomb interactions (a = /3 = 1). Unlike in 
the previous examples, we now find a nontrivial gap func- 
tion that conforms to the "scaling-invariant" momentum 
dependence (15). Computing the free energy for A < fi 
we obtain 



F<{A,^,) 



9 

250 



[/<(A,/i)-/<(Ao,M) 
-IJ^cic^-A^ 



(27) 



where 



7 (A A ln(Ao/c.) [5-ln(AoM] 
/<(A,M) = ^y^ dc^c. ^2 + ln^r/u;r ^''^ 

whereas for A > /i the free energy reads 



2ff( 



F>iA,^,) 

■ [/> (A, fi) - /> (^, ^) + /< (^, ;i) - /< ( Ao, m)] ,(29) 



where 



/>(A,m) 



-A3-2A3ln^ 
3 A 



■In 



Ao 
A 



(30) 



The analysis of the above formulas and the corresponding 
plot of Fig. 3 reveals a local maximum that resides at 
Aoe^^**/^ for /i = and gradually moves rightward as fi 
increases. In contrast, the global minimum located at Aq 
first comes down slightly and then rises with increasing 
/i but, otherwise, retains its position up to /ic ~ O.6A0. 

Upon approaching fic the free energy computed with 
the use of the approximate solution (18) shows a behavior 
suggestive of a first order transition. This observation 
should be contrasted with the conjecture that for = 
the corresponding transition is of infinite order ,5 

The above value of /ic differs from that obtained in 
Ref. 4, in which a first order transition was found at 
/i'c « l-2Ao. It should be noted, however, that in their 
calculation of the free energy, the authors of Ref. 4 in- 
troduced the chemical potential as the lower limit in the 
integration over rather than Ep, which assignment be- 
comes inaccurate for fi Aq. 

Lastly, in Fig. 4 we present the free energy for a = 1/2 
and (3 = 1, which choice of the parameter values corre- 
sponds to the physically relevant problem of the Coulomb 
interacting electrons in an infinite stack of graphene 
layersii^. The behavior exhibited by Fig. 4 is still indica- 
tive of a first order transition at /ic ~ 0.7Ao. 

In fact, a more extensive analysis shows that the above 
pattern appears to be common in the entire range of pa- 
rameters /3/3 < a < f3. The maximum gap Aq (Fig. 5) 
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and the critical coupling gc (Fig. 6), plotted as functions 
of the difference /3 — a agree well with the analytical re- 
sults given by Eqns.13 and 14, respectively. The depen- 
dence of the critical value of the chemical potential /ic 
(Fig. 7) upon that parameter turns out to be somewhat 
less pronounced, though. 

Before concluding, we comment on our assumption of a 
spatially uniform order parameter (that allows one to fo- 
cus solely on the dependence of A on the relative momen- 
tum p of particle-hole pairs). Relaxing this simplifying 
assumption enables one to extend the class of solutions 
and search for possible Fulde-Ferrel-Larkin-Ovchinnikov 
(FFLO) states where the order parameter demonstrates 
periodic spatial variations. Although we choose to post- 
pone a further analysis of this more general case un- 
til future work, our preliminary results suggest that the 
FFLO-like states may only provide a viable alternative 
to the uniform ground state close to the critical value of 
the chemical potential, thus opening the possibility of a 
sequence (at least, two) transitions that might occur in 
the vicinity of the aforementioned ^c- 

As far as the possible applications of the above results 
are concerned, we believe that they may prove helpful for 
resolving a number of outstanding issues pertinent to the 
physics of the nodal strongly correlated electron systems. 

For one, the recent advent of spintronics has brought 
about an upsurge of interest in the magnetic properties 
of nonmetallic compounds, including pyrolytic graphite 
and other degenerate semimetals which can potentially 
harbor a latent excitonic instability responsible for the 
emergence of a (possibly weak, yet robust) ferromag- 
netism, in accordance with the scenario put forward in 
Ref. 7. Although the other recent suspects for the role 
of excitonic ferromagnet, hexaborideSfS, have been largely 
disqualified on the basis of the experimentally established 
extrinsic origin of the observed ferromagnetic response;^ 
the reports of a weak ferromagnetism in some graphitic 
samples^^ has not yet been properly explained. 

Whether the magnetism of graphite is associated 
with a genuine many-body phenomenon such as the 
bulk Peierls (charge and/or spin) instability in gapless 
semimetals or it has a more mundane nature related to 
the presence of (non)magnetic bulk and/or surface de- 
fects, there remains a compelling need for a compre- 
hensive analysis of novel, non-BCS-like, pairing scenar- 
ios in nonmetallic systems with long-range interactions. 
In particular, it is important to ascertain the extent to 
which the customary practice (see Refs. 11 and references 
therein) of replacing the unscreened Coulomb interaction 
with the Hubbard-like on-site repulsion (the only moti- 
vation for which seems to be a greater suitability of the 
latter for numerical simulations) would allow one to un- 
derstand the physics of nonmetallic systems. 

Also, the idea of an incipient phase transition remi- 
niscent of the phenomenon of chiral symmetry breaking 



has recently gained an additional popularity after a num- 
ber of (formally identical but physically distinct) adap- 
tations of the relativistic QED3 theory were discussed in 
such contexts as the antiferromagnet-to-spin liquid and 
pseudogap-to-charge/spin density wave transitions in un- 
derdoped cupratesii^ 

Other examples, including a conjectured second pair- 
ing transition in d-wave superconductors^'^ or /-pairing 
in dichalkogenides;^ have been discussed in the context 
of the (pseudo)relativistic 3D Higgs- Yukawa model. De- 
spite the differences in their underlying physics, all of 
the aforementioned problems appear to be amenable to 
a unified description in the framework of Eq.(4), with 
the parameter values chosen as aqED = /3qed 1 or 
Q!HY = 0, /3hy = 1- Yet another formally related pair- 
ing problem is that of color superconductivity in QCD4 
where the integral kernel in the analog of Eq.(4) behaves 
as U{p — q) oc In |p — . 

The solutions similar to Eq.(15) were also obtained in 
the earlier studies of pairing between incoherent quasi- 
particle excitations mediated by spin fluctuations in d- 
and p-wave superconductors»i^ In that context, the en- 
ergy (albeit not momentum) dependent gap function was 
found to satisfy a counterpart of Eq.(4) with a — f3 — 
1/2. 

In summary, we carried out the analysis of the exci- 
tonic pairing instability in nodal fermion systems, char- 
acterized by a power-law DOS and long-range pairwise 
interaction potentials. We found non-BCS solutions of 
the corresponding gap equation at various finite values 
of the chemical potential. We also computed the free 
energy and studied its dependence on the chemical po- 
tential, thus gaining further insight into the nature of the 
transition in question. We expect these results to be of 
interest for a broad range of physical problems involving 
nodal fermions. 

In particular, we observed that the systems with 
shorter-ranged interactions (lower values of a) and/or 
stronger DOS suppression (higher values of /3) tend to 
manifest a stronger propensity towards excitonic pairing. 
The excitonic order emerges at sufficiently strong repul- 
sive couplings and survives up to a finite critical doping 
by excess carriers. The observed differencies in the be- 
havior at zero and finite a should serve as a warning 
against replacing genuinely long-ranged (e.g.. Coulomb) 
interactions with short-ranged ones, largely for the sake 
of a greater ease of the calculations. 
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FIG. 2: Free energy for q = 0, /3 = 1 as a function of chemical potential /i (in units of Ao): (black), 0.9 (green), 1 (red), and 
1.1 (blue). 
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FIG. 3: Free energy for a = /3 = 1 as a function of chemical potential ^ (in units of Aq): (black), 0.4 (green), 0.6 (red), and 
0.8 (blue). 




FIG. 4: Free energy for a = 1/2, /3 = 1 as a function of chemical potential fi (in units of Ao): (black), 0.5 (green), 0.7 (red), 
and 0.8 (blue). 



